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Abstract 

Spins in quantum dots can act as the qubit for quantum computation. In this context we point 
out that spins on neighboring dots will experience an anisotropic form of the exchange coupling, 
called the Dzyaloshinskii-Moriya (DM) interaction, which mixes the spin singlet and triplet states. 
This will have an important effect on both qubit interactions and spin-dependent tunneling. We 
show that the interaction depends strongly on the direction of the external field, which gives an 
unambiguous signature of this effect. We further propose a new experiment using coupled quantum 
dots to detect and characterize the DM interaction. 
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Solid state devices show great promise for scalable quantum information processing. Sev- 
eral well known proposals for quantum computing have been presented, including semicon- 
ducting quantum dots^ and superconducting Josephson junctions.— Quantum dots currently 
enable the confinement and control of electrons on the scale of tens of nanometers, even 
down to the limit of one electron.— Detection techniques allow the measurement of a single 
electron spin.— 

The solid state matrix provides both opportunities and challenges for quantum control 
and decoherence, due to the complex environment of the qubits. In this paper, we focus on 
a prominent issue for many solid state qubit implementations: the spin orbit interaction, 
which couples spin and charge fluctuations. As typical for qubit interactions, the spin orbit 
coupling can be both useful^ and detrimental.^ Here, we consider how the spin orbit 
coupling affects the time evolution of two-qubit interactions in spin-based quantum dot 
qubits. We find that failure to account for spin orbit coupling can lead to serious control 
errors in the quantum computation. 

The main interaction between spin qubits is the exchange coupling, which can be con- 
trolled with electronic gates, by raising or lowering the electrostatic tunnel barrier between 
neighboring quantum dots^ or by varying the relative depth of the wells constituting the 
double dot.ia I gnoring the crystal matrix, the exchange coupling would be of the Heisenberg 
type, with global SU{2) spin symmetry: -ffHeis = J^i ■ 82- Here, Sj are spin operators and 
J is the tunable exchange coupling constant. The presence of spin-orbit interaction intro- 
duces anisotropy into the exchange coupling, with an antisymmetric component known as 
the Dzyaloshinskii-Moriya (DM) interaction . ^^1^^ Kavokin has shown that the DM exchange 
term also occurs in quantum dots, as a consequence of tunneling coupling.— Devitt et al. 
have proposed methods to determine its magnitude.— Several other authors have also studied 
the importance of the DM interaction for quantum dot quantum computing .^1^^1^^1^^1^^ 

For a single pair of dots, we can write the DM interaction as i^DM = a-^r ■ (Si x S2), 
where r is the unit vector joining the two spins. Thus, the presence of a DM interaction 
reduces the spin symmetry to a global U{1) cylindrical symmetry, where only rotations about 
f remain as symmetry operations. Kavokin has considered the magnitude of the coefficient 
a, computing the exchange integral for the two electrons taking into account the admixture 
of spin projections caused by the spin-orbit interaction.— This calculation applies to the 
case when the energy levels of the individual dots are approximately equal — ^ J) 
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and the Heitler-London method is vahd. Here, we consider only this particular case. In 
GaAs, the predominant spin-orbit coupling is of the Dresselhaus type. For quantum dots in 
a 100 A GaAs quantum well, Kavokin finds a ~ 0.1, not a particularly small value. In such 
a case, we expect the DM contribution to the exchange coupling will be readily apparent. 
For silicon dots, the Dresselhaus interaction is not present, and the predominant spin-orbit 
coupling arises from the Rashba interaction. 

Petta et al. have recently performed a set of experiments with coupled spins in a double 
quantum dot system in GaAs that demonstrate control of the exchange coupling.^'' In these 
experiments, a qubit was defined by \S) and |To): the singlet, and one component of the 
triplet states of the two-spin system respectively. Neglecting additional couplings, we would 
expect the Heisenberg term term to split \S) and |To), thus enabling exchange-based qubit 
rotations.— However, inhomogeneous nuclear fields and the DM interaction also mix in |T+) 
and |T_), the two other components of the triplet. The resulting loss of wavefunction prob- 
ability from the qubit subspace constitutes leakage and it can be interpreted as dephasing 
or decoherence. However, the dynamics are actually coherent. It may therefore be possible 
to utilize the DM dynamics in a beneficial way,'^*^^ or to undo them using time-symmetric 
pulse shapesi^ii^^ or spin echo techniques. Here, we investigate in detail the DM dynamics 
of a double quantum dot system, specifically considering the experiments of Petta et al. We 
explore how the DM interaction modifies the usual interpretation of such experiments, and 
we propose further experiments to detect the presence of the DM interaction and to measure 
its magnitude. 

The Hamiltonian for our double-dot system is 

= J Si • S2 + a J r ■ (Si X S2) - g*^lB (Bi • Si + B2 ■ S2) , (1) 

where Bi = Bext + B„_i and B2 = Bcxt + B„^2- Here, B„^j is the semiclassical field that 
is used to approximate the effective nuclear field for coupling of the electron spin to local 
nuclei in dot i,— given by 

Bn. = ^^El^o(r.)pI^ (2) 

-g ^J'B , 

where l'^ is the nuclear spin operator for a nucleus of total spin / at the lattice site /c, vq is 
the volume of a unit cell containing one nuclear spin, A is the hyperfine coupling strength 
and V'o(^fe) single particle envelope function for the orbital state i evaluated at site k. 
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We assume that B„^j has a Gaussian distribution with mean zero, and a typical variance of 
a = 2.3 mT. In the calculations reported below, our results are averaged over the distribution 
of the nuclear fields,— as consistent with experimental the procedure.— For example, 

F = (27ra2)-3/2 / / desmOx 

Jo Jo 

POO 

/ dB^Ble'^'-Z'^^'FiB^). (3) 
Jo 

High dimensional integrals are evaluated numerically using a simple Monte Carlo integration 
code. 

In Fig. 1, we show the appropriately averaged eigenvalues of H, as a function of the 
applied field -Bext- (Note that we use a = 0.5 in this figure. This large value of a is chosen 
only for purposes of illustration. Elsewhere in the paper we use the more physical value 
a = 0.1.) We observe mixing of the unperturbed singlet and triplet states at special fields. 
Near B^xt = 0, there is mixing of the triplet states, primarily due to inhomogeneous nuclear 
fieldsiS. At nonzero fields, there is an additional mixing of the singlet and triplet states, 
which arises from both inhomogeneous nuclear fields and the DM interaction. The mixing 
occurs near the resonance condition g*fiBBex_t ~ corresponding to Bcxt ~ ±0.04 T in 
the figure. From the point of view of experimental detection, a crucial point is that the 
mixing effect is anisotropic. This is seen clearly in Fig. 1(b) where we plot the overlaps 
between the eigenstates of Eq. ((H) and the pure spin singlet. 

We now compute the time evolution of the two coupled spins for several experimental 
situations of interest. The Hamiltonian in Eq. has four eigenstates \iPi.a), with the 
corresponding eigenvalues -E1..4. For an arbitrary initial state given by |\I/(0)) = 
we can compute the probability Ps(t) = |(S'|\E'(t))p that this state will evolve to a spin 
singlet after time t. We consider the initial state |^I^(0)) = \S). The probability that the 
spin system will remain in its singlet state is then given by 
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P,(t) = ^ 10^1^^ + 25^ |aina,fcos[(Ei-E,)tM]. (4) 

i=l i<j 

Leakage can occur due to both the DM interaction and the inhomogeneous nuclear fields. 
When J is exponentially suppressed, leakage is due entirely to the nuclear fields. For non- 
vanishing J, the initial singlet state would remain stationary if not for the nuclear and 
DM mechanisms. Both mechanisms then play a role in leakage. In Fig. 2, we plot Ps{t) 
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FIG. 1: (a) Energy eigenvalues of the spin Hamiltonian, Eq. as a function of i?cxt; for the 
parameters J = 1 /^eV, a = 0.5 and a = 2.3 mT . Here the sohd hne corresponds to Bcxt||r and 
the dashed Une corresponds to Bext -L r- (b) The overlap |(5'|V'i)P of the energy eigenstates with 
the spin singlet as a function of field. Solid and dashed lines have the same meaning as in (a). 

obtained after allowing the system to evolve over a "waiting time" t = Tg. At the resonance 
condition J = g*fiBBcxt, Ps is strongly suppressed compared to smaller and larger fields. A 
similar suppression of Pg is expected in the absence of spin-orbit coupling. However, the DM 
relaxation mechanism exhibits a strong dependence on the orientation of Boxt with respect 
to r, which cannot be explained by nuclear fields. This dependence on field orientation 
provides an important signature of the DM interaction. 

In the experiments of Petta et al, an initial singlet state is prepared with both electrons 
in a single quantum dot. The electrons are subsequently separated into two dots while 
retaining their singlet correlations. A waiting time ensues, consistent with the analysis 
presented above, after which the singlet probability Pg is measured. The data can be fit 
using a semiclassical model,"21 obtaining a dephasing time of about Tg* ~ 10 ns and effective 
static nuclear field -Bnuc = 2.3 mT. This Tg* is an ensemble averaged time for relaxation to 
the asymptotic value. More relevant for quantum information is the short-time behavior, 
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FIG. 2: (a,b) The probability Ps of an initial spin singlet to remain in the singlet state as a 
function of the "waiting time" for J = 1 /ueV, a = 0.1 and a = 2.3 mT. Solid lines correspond 
to the field orientation Bext||r and dashed lines correspond to Bgxt -L r- (a) -Boxt = 100 mT, (b) 
Bcxt = 40 mT (resonance condition), (c) Dephasing time T*, obtained from a fit to a parabolic 
equation as a function of the orientation angle 6 between Bgxt and r. The fitting data correspond 
to the curves -Bgxt = 40 mT in (b). 

characterized by the quantity T*, defined by P,(r) = 1 - {t/T*Y. In Fig. 2(c), T* is 
plotted as a function of 6, the angle between Bgxt and r. For _Bgxt = 40 mT, the results 
show a significant dependence on 9. Note that the dependence of Ps (as opposed to an 
exponential decay) is due to an absence of dissipation in our model. 

The significance of the DM interaction becomes most apparent during exchange gate op- 
erations, when the Heisenberg and DM couplings, J and a J respectively, are non- vanishing. 
We consider the "Rabi oscillation" experiment of Petta et ai, in which the spins are initially 
prepared in the state |?t,), corresponding to the ground state determined by the nuclear fields 
when J = 0. The initial state is not an eigenstate of Hueis, so when the Heisenberg interac- 
tion is initiated, coherent oscillations will occur between the singlet and triplet manifolds. 
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Thus, after an exchange period of te = 2Tih/ J, the spins will return to their initial state. 
Both inhomogeneous nuclear fields and the DM interaction affect this picture by mixing in 
the different triplet states inhomogeneously, causing P^(t) = \{n\U{t)\n)\^ to decay. Here, 
Pn is the probability to return to the initial state and U{t) is the unitary evolution 

operator for the spin Hamiltonian. If we define the Oj coefficients of the initial state as 
1^) = Yl o.ili'i), then P„(t) is given by Eq. (g)). 

We have computed P„(t) for experimental parameters consistent with Ref. [lO|. The 
results are shown in Fig. 3. Here, the external field is much larger than the nuclear field, 
so the initial state of the evolution is nearly spin polarized. The exchange coupling is 
then switched suddenly to a value slightly off from the resonant condition J = g* HbB cy± 
for a period te- We plot two cases: with and without the DM interaction (a = 0.1,0, 
respectively). In both cases, the initial state |\I'(0)) ~ \Tjr) is very similar to the ground 
eigenstate. (Recall that the nuclear fields and the DM interactions cause a hybridization of 
the IS") and IT4.) states near their level crossing. But away from the crossing, the eigenstates 
retain their \S) and |T+) character.) Therefore, in the long-time limit te ^ Tnuc, Pn does 
not deviate greatly from 1. Here, T^mc ~ ^/fi'*/UB-Bnuc is the nuclear mixing time. We note 
that the solution including the DM interaction is clearly distinguishable from the a = 
case. This is because the DM couphng enhances the hybridization of \S) and |T+), and thus 
the difference between the initial and final states. 

Another obvious feature in Fig. 3 is the initial rapid oscillations of P„. Since the initial 
spin state is not an eigenstate of the exchange Hamiltonian, it can undergo coherent oscilla- 
tions prior to nuclear mixing. In the figure, the predominant oscillations occur between the 
S'-like and T+-like states, with an approximate energy splitting of g*^BBe^t — J and a cor- 
responding oscillation period of 27ih/{g*fiBBcxt — J)- Note that without any true damping 
mechanisms the curves are subject to Poincare recurrence and will return to 1. 

We now propose an experiment to unambiguously detect the presence of the DM inter- 
action. In Fig. 3, the suppression of P„ was strongly enhanced by DM interactions near 
resonance (2.5 /leV). So we perform the previous experiment in a large external field where 
we can tune the exchange coupling to its resonant condition J = g*fiBBcxt during the ex- 
change evolution. Under these conditions, the hybridization of IS*) and |T+) is maximized, 
so that the initial and final states will be quite different. Consequently, after nuclear mixing, 
Pn approaches 0.5. The reason for choosing Pgxt (and thus J) to be large is that this allows 
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FIG. 3: Coherent oscillations of Pn, corresponding to the "Rabi oscillations" of Ref. at 
B = 100 mT, and a = 2.3 mT. Solid curves include the effects of the DM interaction. The dotted 
curve corresponds to DM interactions turned off (a = 0). (A) J = 2.5 /ueV, a = 0.1, (B) J = 3 /ueV, 
a = 0.1, (C) J = 2 /ieV, a = 0.1, (D) J = 3 /ueV, a = 0. 
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FIG. 4: Proposed experiment to observe DM interactions. The system is prepared in the initial 
state |T+) in a large field, -Bext = 1 T. A strong exchange coupling (J = 25.5 /ieV) is initiated at 
the resonance condition J = (7*/UB-Bext with cr = 2.3 mT, producing fast oscillations dampened by 
nuclear mixing (solid curve). In the absence of the DM interaction, no fast oscillations are observed 
(dashed curve). 

many coherent oscillations to occur before nuclear mixing. 

Some typical results are shown in Fig. 4, with and without the DM interaction. Because 
a large value of J has been used, the hybridization of \S) and |T+) is completely dominated 
by the DM interaction for the case a = 0.1. Rapid oscillations occur between these two 
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states, with an energy splitting given by 2|(T'+|ifDM|'S')| = aJ/\/2, and an oscillation pe- 
riod of nVSh/ aJ. To see the fast oscillations, the exchange couphng should be turned on 
quickly compared to the oscillation period, so that the initial state cannot evolve adiabati- 
cally to the ground state. Similar to Fig. 3, the oscillation envelope is eventually suppressed 
by nuclear mixing. However in the large Bcxt limit, the fast oscillation period is deter- 
mined only by DM interactions, not hyperfine effects. This can be confirmed by plotting 
[(fast oscillation period) x B^^i] vs. B^^i, which should remain a constant. The hyperfine 
effects can also be eliminated by polarizing the nuclear spins or by employing a standard 
Hahn spin echo sequence. 

Two-qubit operations require a very accurate knowledge of the spin-spin interaction, and 
the DM interaction is expected to be about a 10% effect in GaAs. It is therefore very 
important to develop methods to measure it in double quantum dot systems. Because the 
interaction breaks spin rotation invariance it can be detected: its effects depend strongly 
on the direction of the applied field in ways that we have described. By carefully choosing 
external parameters, it is also possible to determine the magnitude of the DM coupling by 
measuring the oscillation period for evolution between the singlet and triplet states. 

We gratefully acknowledge conversations with C. M. Marcus and W. O. Puttika. This 
work was supported by NSA and ARDA under ARO contract number W911NF-04-1-0389 
and by the National Science Foundation through the ITR (DMR-0325634) and EMT (CCF- 
0523675) programs. 
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